O: 



O 



SOLUTIONS CONSTRUCTIONS OF A GENERALIZED SYLVESTER 
PROBLEM AND A GENERALIZED FERMAT-TORRICELLI PROBLEM 

FOR EUCLIDEAN BALLS 

^ ■ Nguyen Mau Namjl] Nguyen Hoangj^ and Nguyen Thai aS 

o 

■ Abstract: The classical Apollonius' problem is to construct circles that are tangent to three given 

circles in a plane. This problem was posed by Apollonius of Perga in his work "Tangencies" . The 
Sylvester problem, which was introduced by the English mathematician J.J. Sylvester, asks for 
the smallest circle that encloses a finite collection of points in the plane. In this paper, we study 
the following generalized version of the Sylvester problem and its connection to the problem of 
Apollonius: given two finite collections of Euclidean balls in R", find the smallest Euclidean ball 
that encloses all of the balls in the first collection and intersects all of the balls in the second 
collection. We also study a generalized version of the Fermat-Torricelli problem stated as follows: 
given two finite collections composed of three Euclidean balls in R", find a point that minimizes the 
t~| ■ sum of the farthest distances to the balls in the first collection and shortest distances to the balls in 

, the second collection. 

Key words. Convex analysis and optimization, generalized differentiation, smallest enclosing circle 
problem, Fermat-Torricelli problem. 
^ . AMS subject classifications. 49J52, 49J53, 90C31. 

(N 

1 Introduction 

cn 

I The celebrated Sylvester problem asks for the smallest circle that encloses a finite number 

I of given points in a plane. This problem was introduced by the English mathematician 

James Joseph Sylvester in [1]. Because of its importance in many different applications, 
the problem has attracted many researchers from different fields. The Sylvester problem 
^ ■ and its version in higher dimensions are now called under different names such as: the 

^ . smallest enclosing ball problem, the minimum ball problem, or the bomb problem. The 

readers are referred to [2-7] and the references therein for recent study on the problem and 
its generalizations. 

In the 17th century, the French mathematician Pierre de Fermat introduced another 
problem of Euclidean geometry asking for a point that minimizes the sum of the distances to 
three given points in the plane. The problem was then solved by the Italian mathematician 
and physicist Evangelista Torricelli, and it is now known as the Fermat-Torricelli problem. 
Numerous articles have been written to study this problem and its generalizations to many 
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different settings; see, e.g., [8-12]. 

In our recent publications [13-16], we study extensions of the Sylvester and Fermat- 
Torricelli problems in which the given points are replaced by given sets. The existence 
and uniqueness of optimal solutions, optimality conditions, numerical algorithms, as well 
as other properties of the problems have been initially addressed. In particular, in [16], we 
study of the following problem called the generalized Sylvester problem: given two finite 
collections of sets in a normed space, find a ball of the smallest radius, whose center lies in 
a given constraint set, that encloses all the sets in the first collection and intersects all the 
sets in the second one. We also introduce and study the following generalized version of the 
classical Fermat-Torricelli problem called the generalized Fermat-Torricelli problem: given 
two finite collections of sets in a normed space, find a point in a given constraint set that 
minimizes the sum of the farthest distances to the sets in the first collection and shortest 
distances/distances to the sets in the second collection. 

This paper is a continuation of our development with the study of a special case of the 
unconstrained generalized Sylvester and Fermat-Torricelli problems in which the given sets 
are Euclidean balls. The special features of the Euclidean balls therein make the problems 
distinct from the general case. Although the new problems seem to be very interesting, 
to the best of our knowledge, they have not been considered in the literature. We will 
develop an approach in which the tools of modern convex analysis and optimization are 
employed to solve the problem. Our paper is organized as follows. In section 2, we focus on 
theoretical study of the generalized Sylvester problem for Euclidean balls. An important 
specification for the case of three balls and its connection to the problem of Apollonius are 
carefully investigated. The construction of solutions is presented for each specific problem. 
In section 3, we study the generalized Fermat-Torricelli problem for three Euclidean balls. 
We establish a necessary and sufficient condition for the existence and uniqueness solution 
of the problem in specific cases. The construction of the solutions is also given. 

2 Preliminaries 

In this section, we will formulate mathematical models for the problems under consideration 
and present some concepts and results of convex analysis that will be used in the next 
sections. 

Let / and J be two finite index sets such that |/| + \J\ > 1. Let JB be the unit ball 
and let $7^ := ]B[ai]ri) [ri > 0) for i ^ I and 0j := ]B[bj]Sj) {sj > 0) for j € J be two 
collections of closed balls in R" with the Euclidean norm || • |[. These are our standing 
assumptions throughout the paper unless otherwise stated. 

For a closed, bounded and convex set Q, the farthest distance function and the distance 
function to Q are given respectively by 

M{x; Q) := sup{||x - g|| | g G Q} 

and 

D{x\ Q) := mf{\\x - q\\ \ q e Q}. 
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The generalized Sylvester problem for these Euclidean balls can be reduced to the following 
optimization problem: 

minimize g{x) := max{A^(x), x G M", (2.1) 

where 

M{x) := max{M(2;; r^j) | i € /} and ^(x) := msix{D{x;Qj) \ j € J}. 

Similarly, the mathematical optimization modeling of the generalized Fermat-Torricelli 
problem is 

minimize n{x) := Mi{x) + T>i{x), x G M", (2.2) 

where 

Mi{x) := ^M(x;J7i) and Vi{x) := ^L>(x;ej). 

i&i jeJ 

Let us recall in what follows some important concepts and results from convex analysis 
that will be used throughout the paper. The readers are referred to the books \17\ [18] for 
more systematic development of the field. 

A function : R" ^ R is called convex iff for any x,y G R" and for any a g]0, 1[, one 

has 

'i/'(ax + (1 — a)y) < aip{x) + (1 — a)^(y). 

If this inequality becomes strict whenever x ^ y, the function is called strictly convex. 
A vector v G R" is called a subgradient of a convex function -0 at x iff 

{v, x — x) < tp{x) — iplx) for all x G R". 

The collection of all subgradients of at x is called the subdifferential of the function 
at this point and is denoted by dip{x). In the case where ^/J is Prechet differentiable at x, 
the subdifferential d'il:{x) reduces to the gradient VV'(x) of the function at x, and it is a set 
in the general case. 

The following generalization of the classical Fermat's rule called Fermat subdifferential 
rule will be important in the sequel: 

X is an absolute minimum of a convex function -0 if and only if G dij){x). (2-3) 

Since the functions Q and T-L in problems (j2.ip and (j2.2p are represented respectively in 
terms of the "max" and "sum" of a finite number of convex functions, we are going to use 
available subdifferential rules from convex analysis to further explore these problems. If 
ip{x) := max{'0j(x) | i = 1, . . . , k}, where ipi : R*^ — >■ R for i = 1, . . . , /c are convex functions, 
then 

d^{x) = CO {dtpi{x) I i G I(x)}, (2.4) 

where I{x) := {i = 1, . . . ,k \ ipi{x) = tpix)} is the active index set at x. Another important 
subdifferential rule called the subdifferential sum rule for ijj{x) := X^^Lj^ is stated as: 

k 

d^{x) = ^d^pi{x). (2.5) 

i=l 
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Throughout the paper, we use the following standard notations: bd 17 and int 0, denote 
respectively the boundary and the interior of a set Q; for a,b & M" and a ^ b, 



[a,b] := {ta+{l-t)b [ t G [0, 1]}; 
(a, 6) := (1 - t)b \ t e ]0, 1[}; 

L{a,b) := + (1 - t)b \ t e M}. 



3 A Generalized Sylvester Problem for Euclidean Balls 

3.1 Existence and Uniqueness of Optimal Solutions and Optimality Con- 
ditions 

Let us start the section with simple formulas for computing distances to Euclidean balls in 
M" as well as their subdifferentials in the sense of convex analysis; see, e.g., |17] . 



Proposition 3.1 Let il. = JB{c; r), where c € M" and r > 0. For any x G M", one has 

D{x;n) 



0, if \\x — c\\ < r, 

\\x — c\\ — r, otherwise, 



and 



Moreover, 



M{x;Q) = \\x - c\\ +r. 



N{x;n)ri]B, if \\x - c\\ < r, 
dD{x;n) = { x-c 



otherwise. 



\x 



and 



IB, if x = c, 

dM{x; n) = { x-c 



otherwise. 



X 



For any x £ R", define 

K{x) ■.= {i£l\ M{x;ni) = g{x)} and L(x) := {j £ J \ D{x;Qj) = g{x)}. 

The active index set A{x) is given by the disjoint union A{x) = K[x) U L[x). It is obvious 
that for every x € M", one has 1 < |j4(x)| < |/| + | J|. 

The following theorem establishes a necessary and sufficient condition for the uniqueness 
of an optimal solution. We will provide a simple direct proof for the result; see also [16j . 



Theorem 3.1 The optimization problem (j2.ip always has an optimal solution. Moreover, 
the solution is unique if and only if I ^ ^, or I = % and Dj^jQj contains no more than one 
point. 
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Proof. The fact that problem (j2.ip always has an optimal solution follows from the conti- 
nuity of the function Q therein and the boundedness of its level sets. 

Let us prove the sufficient condition for the uniqueness of an optimal solution. Consider 
the case where 1^0. We will first show that 

0{x) = max{||x — Oill + r^, \\x — bj\\ — Sj : i & I,j G J}. (3.6) 

Since M{x;Qi) = \\x — ai\\ + ri and D{x;Qj) = max{||x — bj\\ — Sj,0} > \\x — bj\\ — sj, one 
always has 

0{x) > max{||x — aj|l + r^, ||x — bj\\ — sj \ i & I,j G J}. 
Fix iQ G /. If K{x) ^ 0, then for an element i € K{x), one has 

Q{x) = \\x — ai\\ + r j < maxjllx — aj|| + ri, \\x — bj\\ — Sj \ i & I, j G J}. 

In the case where K{x) = 0, one has L{x) ^ 0. Fix j € L{x). Then Q{x) = D{x;Qj) > 
M{x; rijo) > 0. Thus x ^ Qj, and hence 

Q{x) = \\x — bjW — rj < maxjllx — ai\\ +ri, \\x — bj\\ — Sj \ i G I,j € J}. 

We have justified (|3.6p . Choose a constant £ such that £ > maxjsj | j G J}. It follows from 
(j3.6p that X is an optimal solution of problem ()2.ip if and only if, it is a solution of the 
following optimization problem: 

minimize G{x), x € M", (3.7) 

where 

g{x) := max{{\\x-ai\\+ri + £f,{\\x-bj\\+£-Sjf \ i G /,j G J} = {g{x)+£f. 

Since fi{x) := {\\x — aj|| + + £)'^ and gj{x) := {\\x — bj\\ + £ — Sj)"^ for i G / and j G J 
are strictly convex functions, we see that problem (|3.7p has a unique optimal solution, and 
hence (12. ip also has a unique optimal solution. 

Suppose that / = and PljgjQj contains at most one point. If Hj^jOj contains exactly 
one point xq, then G^xq) = and xq is the unique solution. In the case Hj^jQj = it 
is also not hard to show that (|3.6p is satisfied for every x G M"", and problem (j2.ip has a 
unique solution. 

Let us now prove the necessary condition for the unique of solution. Suppose that 
problem (j2.ip has a unique solution and assume by contradiction that / = and Dj^jQj 
contains more than one points. It is obvious that g{x) =0 for any x G Hj^jOj. Thus, any 
X G Hj^jQj is an optimal solution of the problem. So we have arrived at a contradiction. 
The proof is now complete. □ 

For any point x G M", the farthest projection and the shortest projection from x to a 
set Q is given by 

V{x; Q) —{q^Ql \\x - q\\ = M{x; Q)} and U{x] Q) := {q £ Q \ \\x - q\\ = D{x; Q)}. 
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If g = iB(c;r), then 



and 

if llx — c|| < r, 

^ix;Q) ; ■ ; ^ ^ ^, . 

■j rr, otnerwise. 

\x — c\\ 

In the case rj = 0, the bah lB{ai;ri) = {aj}, so a ball covers lB{ai;ri) if and only if it 
intersects the ball. Moreover, the case where I = ^ has been considered in [Ij]. Thus, we 
exclude these cases in the theorem below for simplicity. 

Theorem 3.2 Suppose that / 7^ with > for every i £ I. An element x is the optimal 
solution of problem (|2.ip if and only if one of the following conditions holds: 

(1) ]B{x; r), r = Q{x), coincides with k balls from fi, for i £ I , k > 1, lB{x; r) contains the 
other balls in {Qj : i G I}, and it intersects Qj for j £ J. 

(2) V{x;Qi) and Il{x;Qj), i G K(x) and j G L{x), are singletons. Moreover, for pi := 
V{x;Q-i), Qi := JI{x;@j), i G K{x), j G L{x), one has 

X e CO {pi,qi I i G K{x),j G L{x)}. 

Proof. Notice that in this case, problem ()2.ip has a unique solution. By the subdifferential 
Fermat rule, x is the optimal solution of the problem if and only if 

G dg{x) = CO {dM{x;ai),dD{x;ej) \ i G K{x),j G L(x)}. (3.8) 

Let J := {i £ K{x) \ dM{x; Qi) is not a singleton}. If 7^ 0, then x = ai and Q{x) = 
M{x; Qi) = ri for every i £ J'. So all balls Qi for i G ^7 coincide. Moreover, 

ri = M{x; Qi) > M{x; Qj) = \\ai — aj\\ + rj iov i £ J and j £ I \J', 

and Tj = M(x; ilj) > D{x; Qj) = \\ai — bj\\ — Sj for every i £ J and j G J. Thus, ]B{x\ r) 
contains balls in | i G / \ and intersects Qj for j G J. 

Consider the case where J = %. Then dM{x]^i) is a singleton for every i G 
This implies V{x]^i) is a singleton for such an i. For every j G L{x), since -D(x; 0j) > 
M(x;r2jg) > for a fixed index ^ -^j one has that x ^ 0j, and hence JI{x;Qj) is a 
singleton. Then (j3.8p can be equivalently written as 

G CO { -, \ i G K{x),j G L{x)}, where r = = \\x — Pi\\ = \\x — qj\\. 

Thus, there exist Aj > for i G K{x) and fij > for j G L{x) such that X]jg7^(s) Aj + 
T.jem H = 1 and 

n \ " . X — Pi \ ^ X — (Zj 

Equivalently, 

X = ^ AiPi + ^ ^j-gj G CO {pi, qi\i £ K{x),j G L{x)}. 
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The proof of the converse fohows from the Fermat subdifferential rule (|2.3p since we can 
verify that (j3.8p is satisfied in each case. □ 

Finally, using the well-known Caratheodory theorem (see, e.g., |17[ Theorem 1.3.6]) and 
Theorem 13.21 we can prove the proposition below. 

Proposition 3.2 Suppose that I ^ $ with > for every i ^ I . Ifx is the optimal solution 
of problem (|2.ip . then there exist subindex sets h ^ I and Ji C J with 1 < |/i| + |Ji| < n + 1 
such that X is a solution of the generalized Sylvester problem with target sets ili, i G Ii, and 
Oj, j £Ji. 

3.2 Three-Ball Generalized Sylvester Problem and the Problem of 
Apollonius 

In this subsection, we focus on the generalized Sylvester problem for the case of three 
Euclidean balls in two dimensions. From Proposition 13.21 we see that this is one of the 
most important cases since it is possible to reduce the problem with large number of balls 
to the problem of three balls or less. This observation has been used as the key point for 
many algorithms to solve the classical minimum ball problem; see, e.g., [7]. 

Three-Ball Problem: Model I. The first model we study in this subsection is stated 
as follows: given three arbitrary balls = lB{ai;ri) for i = 1,2,3 in the Euclidean plane, 
construct the smallest ball that covers all of the given balls. In this case, / = {1,2,3}, 
J = 0, and problem (j2.ip reduces to 

minimize g{x) = max{M(x; J^i), M{x; O2), M(x; Qs)}, x G (3.9) 

For two balls lB(a; r) and iB(6; s), we say that ]B{a; r) strictly contains ]B(b; s) if IBiJj; s) C 
]B{a;r) and they have no boundary point in common, and ]B(a;r) tangentially contains 
lB{b; s) if lB(b; s) C JB(a; r) and they have exactly one boundary point in common. 

Proposition 3.3 For x gM? and r := G{x), one has that 1^4(3;) | = 1 and x is the solution 
of problem (13. 9p if and only if ]B{x;r) coincides with one of the ball IB{ai;ri) fori G I, and 
lB{x;r) strictly contains the other balls. 

Proof. Observe that problem (j3.9p has a unique solution by Theorem l3.11 Suppose |j4(x)| = 
1, say A{x) = {1}. Since x is the unique solution of the optimization problem (|3.9p . one 
has 

G dM{x;ni). 

Moreover, r = g{x) = M(x;rii), and r > M{x;ni) for i = 2,3. Since G dM{x;ni), 
it follows from the subdifferential representation for M(-;ili) from Proposition 13.11 that 
X = ai, and hence 

r = M{x; Qi) = M(ai; r^i) = n. 

Moreover, r > M(ai,r2j) = \\ai — aj|| + for i = 2,3. This implies lB{x;r) = lB{ai,ri) 
strictly contains two other balls. 
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The proof of the converse is straightforward. Indeed, assume that JB(x; r) coincides with 
]B{ai;ri) and strictly contains the other balls. Then x = ai, Q{x) = r = ri = M(x;J7i), 
and Q{x) > \\ai - ai\\ + u = M{x; Qi) for i = 2, 3. Thus A{x) = {!}, and G dg{x) = B. 
Therefore, x is the solution of problem ()3.9p . □ 

We are going to use the following notations when they are well-defined: 
Ml := P(a2; ^3),vi ■= ^(«3; ^2),wi := ^ ^^ ^i •= r2i); 
M2 := P(a3;ili),U2 := ^(ai; ria), m;2 := ^ , ^2 := 7'(M;2;r22); 
M3 := P(ai;il2),?^3 := P(a2; f^i), -i^s := 1^ , ^3 := ^(if3;^^3)- 




Figure 1: = 2 Figure 2: |^(x)| = 3 

For a,b,c € M^, the notation bac denotes the angle formed by vectors and ah. 

Proposition 3.4 For x G and r := G{x), one has that \A{x)\ = 2 and x is the solution 
of problem (j3.9p if and only if one of the following conditions holds: 

(1) ]B{x;r) coincides with two balls among IB{ai]ri) for i £ I and strictly contains the 
remaining ball; 

(2) ]B{x;r) coincides with one of the balls among lB{ai;ri) for i £ I, ]B{x;r) strictly 
contains another ball, and tangentially contains the remaining ball; 

(3) There exists i £ I such that UitiVi > 90°, x = Wi, and r = ^ ^ — — . 

Proof. Suppose that |j4(x)| = 2, for instance A[x) = {1,2}. Consider the following cases: 
Case 1: If both V{x] and V{x] ^2) are not singletons, then x = ai = 02- Moreover, 

r = M{x\Q.i) = M{x;n2) = n = r2 > M{x;ns), 

so ]B{x;r) = Qi = il2i and lB{x;r) strictly contains Q3. In this case, (1) holds. 
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Case 2: If only one of the sets V{x; Cli),V{x; Q2), for instance V{x] ili), is not a singleton, 
then one has 

X = ai, \\ai — 02!! = n — i'2 and ||ai — asH < n — r^. 

Thus, ]B{x;r) = ]B{ai;ri), lB{x;r) tangentially contains lB{a2',r2), and ]B{x;r) strictly 
contains lB{a^;r^). In this case, (2) holds. 

Case 3: If both pi := 7^(2;; ili) and p2 ■= 'P{x;Q2) are singletons, one has by Theorem 
13.21 that X = 2^^^ • Moreover, ||x — + < _EL__E^_ this case, using the 
representation for farthest projections, we have that x belongs to the open line segment 
(01,02) that connects oi and 02, and x = ws. We also have that pi = us, p2 = V3, and 
ustsvs is greater than 90°. In this case, (3) holds. 

The converse under (1) or (2) follows directly from Theorem 13.21 Suppose (3) holds. 
Then 

M{x;ni) = Af(x;J^2) > Mix^Qs). 

So A{x) = {1,2}. Since 7'(x;ili) = u^, V{x;Q,2) = V3, and x = ^ ^ — by Theorem 13. 2^ 
the element x is the solution of problem (j3.9p . □ 

Proposition 3.5 For x €M? and r := Q{x), one has that \A{x)\ = 3 and x is the solution 
of problem (|3.9p if and only if one of the following conditions holds: 

(1) lB{x; r) coincides with one of the balls among for i ^ I and tangentially contains two 
other balls. 

(2) ]B{x;r) coincides with two of the halls among ilj for i G I and tangentially contains 
the remaining one. 

(3) ]B{x;r) coincides with all three balls. 

(4) V{x;Q,i) are singletons for i G /, x E co {pi,P2,P3}, where pi := V{x;Cli), and r = 
- Pill = P - P2II = 11^ - Pall- 

Proof. Let us prove the implication "if". We have that 

G dg{x) = CO {dM{x; ni),dM{x; J^a), aM(x; ^3)}, 

and r = M(x; J^i) = M{x; Q2) = M{x; ^3). Consider the following cases: 

Case 1: At least one of the sets among V{x;^i) for i € / is not a singleton. If exactly 
one of the sets among V{x] Q-i) for z G / is not a singleton, say V{x] fii), then x = ai and 
r = ri. Moreover, 

M(ai; ilj) = r = |jai — aj|| + for i = 2, 3. 

In this case, Vti tangentially contains two other balls. In the case exactly two sets among 
V{x;Qi) for i ^ I are not singletons, then (2) holds, and if all V{x]^i) for i ^ I are not 
singletons, then (3) holds. 

Case 2: All V{x]Q.i) for i G / are singletons. In this case, one has ||x — Pi\\ = r for z G / 
and by Theorem 13.21 x G co {pi,P2-,P3}- 

The converse follows directly from Theorem 13.21 □ 

We can construct the smallest enclosing ball of three given balls in the plane as follows: 
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Step 1. If there exists one of the three given balls that contains two remaining ones, for 
instance, C f)]^, then x = oi is the solution of the problem and r = ri is the optimal 

value. Otherwise, go to next step. 

Step 2. If one of three angles itjtjUj, i G /, is greater than 90°, then Wi is the optimal 
solution and r = — ^-^ — — is the optimal value. Otherwise, go to next step. 

Step 3. The smallest enclosing ball coincides with the Apollonius ball that is internally 
tangent to ilj for i E /; see, e.g., [19j and the references therein. 

Three-Ball Problem: Model II. The second model we consider in this subsection is: 
given three balls in which are $7j = lB{ai;ri) for i = 1,2 and ©i = lB[bi;si), find the 
smallest ball that covers ili and and intersects ©i. In this case, / = {1,2}, J = {1}, 
and problem ()2.ip reduces to 

minimize g{x) = max{M(x; Jli), M(x; O2), Z)(x; ©1)}, x G M^. (3.10) 

For any u G M^^ one has \K{u)\ G {0,1,2}, [L(n)| G {0,1}, and 1 < \K{u) \ + \L{u)\ < 3. In 
this case, problem ()3.10p has a unique optimal solution by Theorem 13.11 

We say that two balls strictly intersect if they intersect at more than one points, and 
tangentially intersect if they intersect each other at exactly one point. 

Proposition 3.6 For x G and r := g{x), one has that x is the solution of the problem 
(l3T0]l and 

\K{x)\ + \L{x)\ = 1 if and only iflB{x;r) coincides with one of the sets fori=l,2, strictly 
contains the other, and strictly intersects ©1. 

Proof. Suppose that \K{x)\ + \L{x)\ = 1. Then \K{x)\ = 1 and |L(x)| = 0, or \K{x)\ = 
and [^(x)| = 1. 

Case 1: \K{x)\ = 1 and |L(x)| = 0. Suppose that K{x) = {1}. Then 

G dg{x) = dM{x;Qi). 
This implies x = oi, g{x) = M(x; ^li) = ri > M(x; = ||ai — a2|| + 1'2, and 

g{x) = M{x;ni) = ri> Z)(x;0i) = ||ai - - si. 
In this case, we get the conclusion. 

Case 2: |^(x)[ = and \L{x)\ = 1. In this case, one has (7(x) = D{x;@i) > and 

G dG{x) =dD{x;Qi). 
X — bi 

Then x ^ ©i, and hence dD(x; ©1) = -— — -. We arrive at a contradiction since x ^ 61. 

If ~ 

The converse is straightforward. □ 
For {i,j} = {1, 2}, we will use the following notations when they are defined: 

ui := V{bi;ni),vi := n(ai; ©1), = ^"^ '^1 •= ^(^15^2); 

U2 := 'P{bi;n2),V2 := n(a2;©i),w;2 = 1" , ^2 := 'P{w2;^i)- 
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Proposition 3.7 For x € and r := Q{x), one has that x is the solution of the problem 
(j3.10p and \K{x)\ = \L(x)\ = 1 if and only if one of the following conditions hold: 

(1) ]B{x;r) coincides with one of the sets 0,i for i=l,2, strictly contains the other, and 
tangentially intersects Qi. 

(2) One of the angles UitiVi for i = 1,2 is greater than 90°, x = Wi, and r = ■ 



2 

Proof. Suppose that x be a solution of the problem and |-R'(x)| = |-L(x)| = 1. Then 
K(x) = {1} or K{x) = {2} and L{x) = {1}. Consider the case where K(x) = {1} and 
L(x) = {1}. Then 

r = g{x) = M{x; Oi) = D{x; Gi) > M{x; fis) 

and 

G dg{x) = CO {dM{x; Qi),dD{x; Qi)}. 

Since D{x; Qi) > 0, one has x ^ ©i. Observe that if P(x; is not a singleton, then x = ai, 
and hence ri > M(x; 0,2) = Wo-i — a2\\+r2, and ri = M{x; Oi) = D{x;@i) = ||ai — 61 1| — Si. 
In this case, (1) holds. 

Suppose that 7^(x;0i) be a singleton. By Theorem 13.21 x = — ~; where yi = 
V{x; and 

zi = n(x;©i). In this case, yi = ui, zi = vi, and x = wi. Since r = ^ — — > 

M{wi;0,2) = \\wi — till, one has that uitiVi > 90°. The proof of the converse is also 
straightforward. □ 

Denote 

Xl ~\~ X2 

Xl := V{a2;0.i),X2 := V{ai;Q2),y = ,z := U{y;&i). 



Proposition 3.8 For x € and r := G{x), one has that x is the solution of the problem 
(jS.lOp and \K(x)\ = 2, \L{x)\ = if and only if one of the following conditions hold: 

(1) ]B{x;r) coincides with both for i = 1,2, and strictly intersects Qi. 

(2) ]B{x;r) coincides with one of the balls fij fori = 1,2, tangentially contains the other, 
and strictly intersects ©i. 

(3) The angle X1ZX2 is greater than 90°, x = y, and r = ^- — 

Proof. Suppose that x be a solution of problem (j3.10p . K{x) = {1, 2}, and L(x) = 0. Then 
r = g{x) = M{x,ni) = M{x,Vl2) > D{x,Qi), 

and 

G dg{x) = CO {dM{x; ^i), dM{x; Q2)}- 
Let us consider the following cases: 

Case 1: If x = ai and x = 02, then ri = r2 = r, and hence IB{x;r) = Qi = 0,2- On the 
other hand, since D[x, ©1) < r , one has that ]B(x;r) strictly intersects ©1. 
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Figure 3: \K{x)\ = 1, \L{x)\ = 1 Figure 4: \K{x)\ = 2, \L{x)\ = 

Case 2: If x = ai and x ^ a2, then ri = r = r2 + — a2||. In this case, we have 
]B{x;r) = Since — a2|| = — r2 and D{x,@i) < r, we have that lB{x;r) tangentiahy 
contains and strictly intersects ©i. 

Case 3: Now, we consider the case x ^ ai and x ^ a^- Using Theorem 3.2, we see that 

X = txi + (1 - t)x2, t G ]0, 1[. 

X\ ~t~ X2 

Since the problem has a unique solution, then t = 1/2, i.e., x = ^ = y and r = 

\\xi — X2\ 



2 . Since D{x, ©i) < r, or Hy — z|| < r, the angle 2:1^x2 is greater than 90°. 
The proof of the sufficient condition is straightforward. □ 

Proposition 3.9 For x € and r := G{x), one has that x is the solution of the problem 
()3.10p and \K(x)\ = 2, \L{x)\ = 1 if and only if one of the following conditions hold: 

(1) ]B{x\r) coincides with both of the halls Vti fori = 1,2, and tangentially intersects ©i. 

(2) lB{x;r) coincides with one of the halls Vti fori = 1,2, tangentially contains the other, 
and tangentially intersects ©1. 

(3) V{x; rij) for i = 1,2 and n(x; ©1) are singletons, \\x —pi\\ = ||x — P2II = ||^ — where 
Pi = V{x; Qi) for i = 1,2 and qi := Il{x; ©1), and x G co {pi,P2, Qi}- 

Proof. Suppose that x be a solution of the problem and K{x) = {1, 2} and L{x) = {1}. 
Then 

r = g{x) = M(x, Oi) = M(x, ^2) = D{x, ©1), 

and 

G dg{x) = CO {dM{x; Qi),dM{x; Q2),dD{x, ©1)}. 
We consider the following cases: 

Case 1: If x = oi and x = 02, then ri = r2 = r, and hence JB{x;r) = Qi = Q.2. On the 

other hand, since D[x, ©1) = r, one has that ]B{x;r) tangentially intersects ©1. 

Case 2: If x = ai and x 7^ 02, then ri = r = r2 + \\x — a2\\. In this case, we have 
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lB{x; r) = ill. Since ||x — a2|| = n — r2 and Gi) = r, one has that r) tangentiahy 
contains 0.2 and tangentially intersects 0i. 

Case 3: Now, we consider the case where x ^ a\ and x ^ a2. Then for i = 1,2 

are singletons. Since D(x,0i) = r > 0, n(x,Bi) is also a singleton. 

Using Theorem 13.21 we see that x = tipi + t2P2 + ^3^1, ti G [0, 1], where ti + ^2 + ^3 = 1- 
Moreover, ||x — pi || = ||x — P2II = ||^ ~ 91 II- 

The proof of the sufficient condition is straightforward. □ 




Figure 5: \K{x)\ = 2, |L(x)| = 1 

We are now able to construct the smallest ball that corresponds to the solution of 
problem (j3.10p as follows: 

Step 1. If there is a ball from ^li and ^12 that contains the other and intersects ©i, then 
that ball is the solution of the generalized Sylvester problem. Otherwise, we go to next 
step. 

Step 2. One of the angles uttiVi for z = 1, 2 is greater than 90°, x = Wi, and r = ' ^ — —. 
Otherwise, we go to next step. 

\\x\ — -^2 

Step 3. The angle xizx2 is greater than 90°, x = y, and r = . Otherwise, we go 

to next step. 

Step 4. In this case, the smallest ball is the Apollonius ball that is internally tangent to 
0,1,^12 and externally tangent to O. 

Three-Ball Problem: Model III. The third model we will consider is: given three balls 
0,1 = ]B{ai;ri), ©i = lB(bi;si), and ©2 = IB(b2',S2), find the smallest ball that covers Oi 
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and intersects 0i and ©2. In this case, / = {1}, J = {1,2}, and problem (j2.ip reduces to: 

minimize g{x) = max{M(x; ni),D{x; ©i), D{x; G2)}, x G M^. (3.11) 

Proposition 3.10 For x G and r := g{x), one has that x is the solution of the problem 
()3.1ip and \K{x)\ + \L{x) \ = 1 if and only if ]B{x;r) coincides with $7i and strictly intersects 
both ©1 and ©2. 

Proof. Suppose that \K{x)\ + \L{x)\ = 1. Then \K{x)\ = 1 and |L(x)| = 0, or \K{x)\ = 
and = 1. Consider the first case where = 1 and = 0. Then K{x) = {1} 

and 

G dM{x]ni). 

Moreover, r = M{x; Vli) > D{x; ©j) for i = 1,2. This implies x = ai, and hence 

M{x; Qi) = ri = r> D{x; ©i). 

In this case, lB{x;r) coincides with Qi and strictly intersects ©1 and ©2- 

For the case where = and = 1, we can assume that L{x) = {1}. Then 

G dD{x;@i), 

D{x;Qi) > Z)(x;©2), and D{x;Qi) > M{x;Ui). This implies i:>(x;©i) > 0, and hence 
X ^ ©1. We have arrived at a contradiction due to the representation of the subdifferential 
of distance function at out-of-set points from Proposition I3.1[ The converse follows from 
Theorem [321 □ 
For {i,j} = {1, 2}, we will use the following notations when they are defined: 

c- + d- 

a = V{hi;Qi),di = n(ai;©i),ei = ' ^ \ fi = II(ei;©j). 



Proposition 3.11 For x G and r := G{x), one has that x is the solution of problem 
(13. lip and \K{x) \ = \L{x) \ = 1 if and only if one the following holds: 

(1) ]B(x;r) coincides with Qi, tangentially intersects one of the ©j fori = 1,2 and strictly 
intersects the other. 

(2) There exists i £ J such that Cifidi > 90° , x = Ci, and r = — ^ — 

Proof. Suppose that x be a solution of the problem and |-f^(x)| = |-^^(x)| = 1. Then 
K{x) = {1} and [L{x) = {1} or L{x) = {2}]. Consider the case where K{x) = {1} and 
L{x) = {1}. Then 

r = g{x) = M{x; Oi) = D{x; ©1) > D{x; ©2) 

and 

G dg{x) = CO {dM{x; Qi),dD{x; ©i)}. 
Since D{x; ©1) > 0, one has x ^ ©1. 
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In the case 7^(3;;Oi) is not a singleton, one has x = oi. Then ri = M{x;Qi) = 
D(x; ©i) = ||ai — 61 II — si and ri = Af(x; ^li) > D{x; ©2). In this case, (1) holds. 

Now assume that T'{x; is a singleton. By Theorem 13.21 one has x = — — where 
yi = V{x;Qi) and zi = n(x;©i). In this case, zi = ci, yi = di, and x = ei. Since 
r = ||ci — (ii||/2 > D{ei;@2) = ||ei — /i||, one has that Cifidi > 90° for i = 1. The converse 
follows from Theorem 13. 2[ □ 




Figure 6: \K{x)\ = 1, \L{x)\ = 1 Figure 7: \K{x)\ = 0, |L(x)| = 2 

For i,j £ {1, 2}, i 7^ j, we will use the following notations when they are defined: 
yj = U{bi, ©,), z = t = V{z; fti). 

Proposition 3.12 For x G and r := G{x), one has that x is the solution of the problem 

_ _ ^ — ^ _ \\)\ — y2|| 

(|3.1ip and \K{x)\ = 0, \L{x)\ = 2 if and only if yity2 > 90°, x = z, and r = . 



Proof. In this case, K{x) = and L{x) = {1,2}. Then 

G dg{x) = CO {dD{x; ©1), dD{x; ©2)}, 



and r = D{x;Qi) = Z)(x;©2) > M{x;0,i). This implies x ^ ©j for i = 1,2. Let Zi 
n(x; ©i) for i = 1 
X = z. Moreover, 



n(x;©j) for i = 1,2. By Theorem 13.21 x = — -. This implies Zj = yi for i = 1,2 and 



m - y2 ^ n ^ II II 

> M[x; = \\ai — z\\. 



2 

Thus, yity2 > 90°. The proof of the converse is straightforward. □ 
The proof of the proposition below is similar to that of Proposition 13.91 

Proposition 3.13 For x € and r := Q{x), one has that x is the solution of the problem 
(|3.1ip and \K{x) \ = 1 and L{x)\ = 2 if and only if one the following holds: 

(1) ]B{x;r) coincides with ili, is tangent (externally) to both ©j for i = 1,2. 
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Figure 8: \K{x)\ = 1, \L{x)\ = 2 



(2) V{x; n(x; Qi) for i = 1,2 are singletons, \\x — pi|| = ||a; — 9i|| = — 92II) where 
pi = V{x; and qi := Il{x; Qi) for i = 1,2, and x G co {pi, qi,q2}- 

We are now able to construct the smallest ball that corresponds to the solution of 
problem (j3.1ip as follows: 

Step 1. If rii n 01 7^ and 17i n 02 7^ then x = ai is the solution of the problem and ri 
is the optimal value. Otherwise, go to next step. 

Step 2. If one of the angles Cifidi for i = 1, 2 is greater than 90°, then x = Ci, r = — 
Otherwise, go to the next step. 

— _ llyi ~ 2/211 

Step 3. If the angle yity2 > 90°, then x = z and r = . Otherwise, go to the next 

step. 

Step 4. In this case, the smallest ball is the Apollonius ball that is internally tangent to 
0,1 and external tangent to 0i,02. 

Three-Ball Problem: Model IV. Let us now consider the smallest intersecting ball 
problem: given three balls 0j = lB{hi; Si) for i = 1,2, 3, find the smallest ball that intersects 
0j for i = 1,2, 3. In this case, / = 0, J = {1, 2, 3}, and problem ()2.ip reduces to 

minimize g{x) = max{Z)(x; Qi),D{x; Oa), D{x; 63)}, x G (3.12) 

In the case n^^j^Gj 7^ 0, any point in this intersection is a solution of problem (j3.12p . so we 
only consider the case where this intersection is empty. By Theorem 13.11 problem (j3.12p 
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has a unique solution. It is also not hard to see that = |^(x)| > 2. 

We will use the following notations when they are defined: 

ui = [62, 63] n bd{Q2), VI = [62, h] n K03), 

U2 = [61,63] n6(i(e3), V2 = [61,63] n6(i(ei), 
^3 = [61,62] n6(i(ei), 1^3 = [61,62] n6(f(e2), 

Ui +Vi U2+ V2 U3 + ^3 

mi = , m2 = , ms = , 

xi=n(mi;6i), X2 = n(m2;02), X3 = n(m3;03), 

= Birrn; ^^1^), B2 = B{m2; ^^^), B, = B{m,- 



/ S 
/ \ 




Figure 9: \A{x)\ = 2 Figure 10: |^(x)| = 3 

It is not hard to see that the smallest intersecting ball can be constructed as below: 

Step 1: If there exists j E {1,2,3} such that uJxjVj is greater than 90°, then Bj is the 
solution of the problem. Otherwise, go to the next step. 

Step 2: The smallest intersecting ball is the Apollonius ball that is externally tangent to 
the three balls. 

The readers are referred to |14) for a similar construction for the case where three given 
balls are disjoint. 

4 Generalized Fermat-Torricelli Problem for Euclidean Balls 

In this section, we consider the following problem: 

minimize n{x) = M{x; a) + D{x; Sj), x G R"^, (4.13) 
iei jeJ 
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where |/| + | J| = 3. It is clear that x is an optimal solution of problem (j4.13p if and only if 
it is a solution of the following problem: 

minimize T-L{x) = ||x — ai|| + D{x; Qj), x G M^. 
Since a singleton is a ball with radius 0, it suffices to consider the following problem: 

3 

minimize U{x) = ^ D{x; Qi), xe M^. (4.14) 

i=l 

Throughout this section, we assume that the given balls have distinct centers since the 
other case is trivial. 

4.1 Existence and Uniqueness of Optimal Solutions 

In this subsection, we will study properties of solutions of problem (j4.14p that enable us to 
derive a necessary and sufficient condition for the problem to have a unique solution. 

Proposition 4.1 The solution set S of problem (j4.14p is a nonempty compact convex set 
in E^. 

Proof. It is clear that the function 1-L{x) is continuous. For any 7 > inf2.g]g2 'H{x), we see 
easily that {x G | < 7} is compact. Thus, (I4.14p always has an optimal solution. In 

particular, S is compact. Since H is convex and continuous, the solution set of the problem 
is also convex. □ 
For any u gM?, define 

A{u) := {i e {1,2,3} | u € 9,}. 

Proposition 4.2 Suppose that \A(x:^)\ = 0. Then is an optimal solution of problem 
(j4.14p if and only if x^, is the solution of the classical Fermat- Torricelli generated by the 
centers of the balls: 61,62,63. 

Proof. Suppose |A(x*)[ = and x* be an optimal solution of problem (I4.14|) . Then x^ ^ Qi 
for i = 1,2, 3. Choose 6 > such that ]B{x^,; (5) H ©j = for i = 1,2, 3. Then the following 
holds for every x € IB{x^; 5): 



3 

' Si 

1=1 i=l ^ i=l i=l 



||x* - 6i|[ - ^ Sj = inf 7i{x) = ?^(x*) < 7i{x) = ^ ||x - 6j|| - ^ 

1=1 i=l ^ 

Thus, X* is a local minimum of the problem 

J^{x) := ||x — 6i||, X € 



3 

mmimize ] 

1=1 



so it is also an absolute minimum of the problem since is a convex function. Therefore, 
X* the solution of the classical Fermat-Torricelli generated by 61 , 62 , 63 . The converse is also 
straightforward. □ 
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Lemma 4.1 Let S he the solution set of problem (|4.14p . Suppose that there exists E 5 
such that A{x^) = 0. Then S is a singleton, namely S = 

Proof. Suppose that A{x^,) = 0, and there exists £ S with x* 7^ y*. Since A{x^:) = 0, 
one has x* ^ Qi for i = 1,2,3. Then x^, is the solution of the classical Fermat-Torricelli 
problem generated by the centers of the balls: 61,621 and 63. Since S is convex, one has 
Q S. Thus, it is possible to find 7^ x*, G S, and z* ^ 0j for f = 1,2,3. Then 
z=K is also the solution of the classical Fermat-Torricelli problem generated by the centers of 
the balls, which is a contradiction since this problem has a unique solution. □ 

Lemma 4.2 For any a > 0, and a, 6 € with a ^h, consider the set 

:= {x e I ||x - o|| + ||x - 6|| = a}. 

X ~\~ y X -\- y 

Suppose that \\a — b\\ < a. For x,y £ E with x ^ y, one has \\ — a\\ + || — 6|| < a, 

X ~\~ y 

and in particular, — - — ^ E. 

Proof. We have 

uX + y 11 ..x + y 1 

11^ all + 11^ 6|| =-(|[x - a + y - a|| + ||x - 6 + y - 6||) 

<-(||x - a|| + ||y - a|| + ||x - 6|| + \\y - 6||) = a. 

Suppose by contradiction that 

11^ all + 11^ 6||=a. (4.15) 

X ~\~ y 

That means — - — G E. By a property of the Euclidean norm, one has 

X — a = k{y — a) and x — 6 = m{y — 6) 

for some numbers k, m in (0, +00) \ {1} since x ^ y. This implies 

1 k 1 m 

-X - ry e L(x, y) and 6 = x - y G L(x, y). 



1 — k I — k 1 — m 1 — m' 

X ~\~ y 

Since a > ||a — 6||, it is not hard to see that x,y E L(a, 6) \ [a, 6]. One also has — ^ — G [a, 6]. 

X ~\~ y 

Indeed, assume for instance, that the order of points is: x, — ^ — , a, 6, y. Then 

II — 2 '^ll " — 2 ~ ~ ^ 

which contradicts ()4.15p . Now, since — ^ — € [a, 6], one has 

uX + y + y 

II — a\\ + II — o|| = ||a — o|| < a. 
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which is a contradiction. We have proved that || 0^11 + II ^11 < ct- Thus, 

^^E. 2 2 

2 ^ 

Let be a subset of M". We say that O is strictly convex if for any x,y € Q with x ^ y 
and for any t S ]0, 1[, one has tx + {1 — t)y G int 0. In the setting of Lemma 14.21 the set 

E:={xeR'^ I llx - a|| + ||x - 6|[ < a} 

is strictly convex. This is in fact a weh-known result, but we provide the detailed proof for 
the convenience of the readers. 

Lemma 4.3 Let S be the solution set of problem (I4.14p . Suppose that there exists S 5" 
such that A{x^:) = {i} and ^ [bj,bk], where i,j,k are distinct indices in {1,2,3}. Then 
S is a singleton, namely S = {x*}. 

Proof. Suppose without loss of generality that A{x^) = {!}. Then and G Qi, x^ ^ ©2, 
X* ^ 63, and X* ^ [b2,bs]. Let 

a := ||x^, — 62II + ||x* — 63II = inf Hix) + S2 + S3. 

Then a > \\b2 — 63 1|. Define 

£' := {x G I ||x - 62II + II2; - 63II = a}. 

Clearly, x* S E H @i. Suppose by contradiction that S be not a singleton. Then there 
exists y* G S and 7^ x*. Since S is convex, [x*,y*] C S. We can choose G that 
is close enough and distinct from x* such that [x^,, z^] n ©2 = and [x*, z*] n ©3 = 0. Let 
us first show that ^ ©1. Indeed, assume that G ©i. Then D(z*;©i) = 0, and hence 

7i(z^,) = inf ?^(x) = D(z^; ©2) + Z?(z*; ©3) = |k* - ^2|| - S2 + |U* - ^sll - S3. 

This implies z^, G ii^. Since S and ©1 is convex, —-^ — - £ S" n ©1. Moreover, * ^ — ^ ^ ©2 

and * ^ — ^ ^ ©3. Again, we have that * ^ — ^ G i?. This is not the case by Lemma |4.2[ 
We have shown that z* G 5" and A{z^) = 0, that contradicts the result from Lemma HTJ 
Therefore, S must be a singleton. □ 

Lemma 4.4 Let S be the solution set of problem (|4.14|) . Suppose that x* G S* and yl(x*) = 
{i,j} C {1,2,3}. Thenx^ G bd (©in©^). 

Proof. Assume without loss of generality that A{xtf) = {1,2}. Then x* G ©1 fl ©2 and 
X* ^ ©3. Suppose by contradiction that x,,, G int (©1 n ©2). Then there exists 5 > with 
iB(x*; (5) C ©1 n ©2. Let p := n(x*; ©3) and 7 := ||x* — p|| > 0. Let q G [x*,|)] fl iB(x*, 5) 
satisfy ||(? — p|| < ||x^,— p|| =D(x^,;©3). Then 

3 3 
n{q) = D{q; e,) = D{q- ©3) < \\q-p\\ < \\x,-p\\ = D{x,; ©3) = D{x,; Qi) = H{x,). 

i=l i=l 

This is a contradiction to the fact that x* G 5. □ 
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Theorem 4.1 Suppose Cif^iQi = 0. Problem (j4.14p has more than one solution if and only 
if there is a set [bi,bj] n G^, for distinct indices i,j,k G {1,2,3}, that contains a point u 
belonging to the interior of Qk CLnd \A{u)\ = 1 . 

Proof. Suppose that [61,^2] contains a point u that belongs to the interior of ©3 and 
[^(u)! = 1. Then n ^ ©1 and u ^ ©2. We can choose f 7^ n, v & [&i,^2]> v ^ ©1, v ^ Q2 
and V G int03. Then |^(v)| = 1. We first prove that u is a solution of the problem. Indeed, 
in this case, 

niu) = D{u- ©i) + D{u] ©2) + D{u; ©3) = - 62II - si - sg. 

Choose (5 > such that B{u]5) n ©1 = and JB(n;5) n ©2 = 0, and B{u;5) C ©3. For 
every x € ]B{u;5), one has 

7i{x) = D{x; ©2) + D{x; ©3) = \\x - bi\\ + \\x - 62II - si - ■52 > ll^i - ^2!! - si - ■52 = Ti{u)- 

This implies that n is a local minimum of 7^, so it is also an absolute minimum since 'H(n) 
is convex. An analogous argument can be applied for v. Then u, v S. 

Now assume that S has more than one elements. Let be two distinct elements 

of S. Then [ 3^*5 2/*] ^ S by Proposition UlTl If there is a solution that does not belong to 
©j for every i G {1,2,3}, then S reduces to a singleton, that is a contradiction. We can 
assume without loss of generality that ©1 contains infinitely many solutions. Then int ©1 
contains infinitely many solutions by the strict convexity of ©1. If there is such a solution 
u with A{u) = {1}, then u G int ©1 and u does not belong to ©2, ©3. Thus, u G [62,^3]) 
since if not, the solution must be unique by Lemma 14.31 Therefore, the conclusion holds. 
Suppose = 2 for every solution that belongs to int ©1. Then there are infinitely 

many solutions that lies on the intersection of two sets, which is strictly convex in this case. 
So there is a solution that belongs to the interior of this intersection, that is a contradiction 
to Lemma |4.4[ □ 

Example 4.1 In Figure 1, the given generalized Fermat-Torricelli problem has infinite 
solutions because [&i,^'2] intersects the interior of ©3, and = 1 for every u G {B,C). 

In fact, the solution set is S = [B,C]. 

Let Vk := [bi,bj] flint ©^ for distinct indices i,j,k G {1,2,3}. As a direct consequence of 
Theorem 14.11 we have the following corollary. 

Corollary 4.1 Suppose nf^^©i = 0. Problem (j4.14p has a unique solution if and only if 
the following implication holds for any k G {1, 2, 3}.' 

[Vk / 0] ^ = 2 for ah n G Vfc] . 

The corollary below provides a sufficient condition for problem ()4.14p to have a unique 
solution. 

Corollary 4.2 Suppose nf^^Oi = 0. Problem (j4.14p has a unique solution if [bi,bj] fl ©^ 
contains at most one point for distinct indices i,j, k G {1, 2, 3}. 

Proof. Suppose by contradiction that (j4.14p has more than one solutions. By the previous 
theorem, there is an interval, say [61,62]) that contains a point belonging to the interior of 
©3. Then [61,62] H ©3 contains infinitely many points, that is a contradiction. □ 
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Figure 11: A generalized Fermat-Torricelli with infinitely many solutions 
4.2 Solution Constructions 

In this subsection, we will propose a method of constructing a solution of problem ()4.14jl 
for arbitrary balls Gj, i G J = {1,2, 3}. 

Proposition 4.3 An element is an optimal solution of problem (j4.14p if and only if 

- Yl ""i^ Yl [N{x.;@i)nB], (4.16) 

i€J\A(x,) i£A{xt) 

where e,- := 



X* bi 



Proof. By the Fermat subdifferential rule (j2.3p . x=k is an optimal solution of problem (j4.14p 
if and only if 

ieJ i€J\A{x,) jeA(i'*) 

The conclusion then follows from Proposition 13.11 □ 
Proposition 4.4 Consider problem ()4.14p . Suppose that nf^^Qi 7^ 0. Then 

s = niiGi. 

Proof. Fix any x* € flf^^Gi 7^ 0. Then 7i{x^) = 0, and hence x^, G 5 since 

inf nix) > 0. 
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Conversely, for any x G 5, one has ^{{x) < T-L{x:^) = 0. Thus, D{x;Qi) = for i = 1,2,3. 
This imphes x G @i for i = 1, 2, 3, or equivalently, x G r\f^i@i 7^ 0. □ 

From now on, we only need to consider the case where n^^^Qj = 0. Then < 3 

for all X:^ € S. 

Proposition 4.5 Suppose that = 2 andx^, be an optimal solution of problem (I4.14p . 

say X* € ©1 n 02- Then one of the following condition holds: 

(1) x^: is the intersection of [61,63] and the boundary o/Gi, and G int 02- 

(2) x^: is the intersection of [62763] CLnd the boundary 0/G2, and x* € int Gi. 

(3) x^, belongs to the intersection of the boundary of Qi and the boundary of Q2, o-'^d 

—63 = sei + ie2, 
where s,t £ [0, 1] and Cj = -r^^ — ^ for i = 1,2, 3. 

1 1 1 1 

Conversely, if one of the conditions above is satisfied, then x* is an optimal solution of 
the problem. 

Proof. Suppose that |^(x*)| = 2 and x* G Gi n ©2. By Proposition 14.3^ 

- 63 G [N{x,; ©i) n iB] + [iV(x,; G2) n B]. (4.17) 

Since 63 ^ 0, one has that x* G bd ©1 or x* G bd ©2- 

Case 1: x* G bd Gi and x* G int ©2. In this case, (|4.17p reduces to 

-63 G iV(x*;©i). 

This is equivalent to the fact that x* is the intersection of [61, 63] and the boundary of ©1. 
Thus, (1) is satisfied. 

Case 2: x^, G bd ©2 and x^, G int ©1. Similar to Case 1, under this condition, (2) is 
satisfied. 

Case 3: x* G bd ©1 and x* G bd ©2- In this case, (j4.17p reduces to 

—63 = S61 + te2, 

where s, t G [0, 1]. Notice that in this the intersection of bd ©1 and bd ©2. Since 

there are at most two point in this intersection, this condition is verifiable. □ 

Proposition 4.6 Suppose that \A{x^)\ = 1 andx^ be an optimal solution of problem (j4.14p . 
say X* G ©1. Then 

-62 - 63 G A^(x*; ©1) and (62, 63) < -1/2. 
Conversely, if this condition is satisfied, then solution of problem (I4.14p . 

Proof. Condition (j4.16p can be written as 

-62 - 63 G A^(x*; ©1) and ||62 + 63II < 1. 
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Notice that ||e2 + e^W < 1 if and only if (62,63) < —1/2. The proof is complete. □ 
The construction to find a solution is as follows. We first will find all possible solutions 
in each set Gj for i = 1,2, 3. If a solution is found, then there is no solution outside of the 
sets by Lemma 14. 1[ If no solution is found, then the optimal solution is found by solving 
the classical Fermat-Torricelli generated by bi for i = 1,2,3. For instance, solutions on Qi 
can be found by the following steps: 

Step 1: If nf^^Qi 7^ 0, then every point in this intersection is a solution. Then, go to the 
next step. 

Step 2: Connect the centers j = 2,3. If, for instance, [61,^3] intersects bd ©i at a 

point that belongs to int ©2, then that point is a solution. Then, go to the next step. 

Step 3: Find the intersection of the boundary of the balls ©1 and ©2, and ©1 and ©3. For 
instance, let p and q be the intersections of bd ©1 and bd ©2. Then verify if the condition 
—63 € [0, l]6i + [0, 1]62 is satisfied at each point. If the condition is satisfied, for instance, 
at p, then p is a solution. Then, go to the next step. 

Step 4: Verify if [62,^3] intersects ©1. If, for instance, [62, ^3] intersects ©1, then find all 
point in the intersection does not belong to ©2 and ©3, and those are solutions. In the case, 
[62, ^3] does not intersect ©1, then we find U2 ■= [bi, 62] fl bd ©1, and U2 ■= [bi, 63] n bd ©1. 
Next, find a unique point qi on the minor curve ^2^3 such that b2qibi = ^39161. If qi does 
not belong to ©2 and ©3, and ^291 ^3 ^ 120°, then qi is a solution of the problem. 
The same process can be repeated to find solutions on ©2 and ©3. 

5 Concluding Remarks 

In this paper, we have provided a detailed theoretical analysis for the generalized Sylvester 
problem and the generalized Fermat-Torricelli problem for Euclidean balls. A natural ques- 
tion is: can we develop numerical algorithms to solve the generalized Sylvester and the 
generalized Fermat-Torricelli problems for Euclidean balls proposed in this paper? The 
well-known subgradient method provides such a simple algorithm since subgradients of the 
distance functions involved in the problems can be explicitly determined. By exploiting the 
majorization-minimization (MM) principle of computational statistics, another algorithm 
for solving the constrained version of the generalized Fermat-Torricelli has been developed 
in [20]. Faster algorithms can also be developed using other methods of nonsmooth opti- 
mization. We will address this in our future research. 
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